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These transfer functions then represent the effects of the 
fundamental solar array bending modes on the rotational 
response of the rigid bus. One additional simplification 
was made before incorporating these models into the ana- 
log computer simulation of vehicle dynamics. Since the 
quantity (1 — 96.64/1*^) is very nearly 1, meaning that the 
structural pole and zero contributed by the solar array to 
the (/-axis loop are very nearly coincident, virtually no 
structural effects on control response should be expected 
from this mode. It was therefore eliminated. 

The complete analog computer simulation program for 
the attitude dynamics (including scaling) is shown in 
Fig. 13. The simulation required the use of five electronic 
(diode “quarter-square”) multipliers, 13 integrators, 36 
summing or inverting amplifiers, approximately 35 poten- 
tiometers, and a Brush Instruments strip-chart recorder. 


F. Some Closed-Loop TVCS Test Responses 

Figures 14 and 15 show the recorded responses of sev- 
eral variables within the TVCS (modified as discussed in 
Sections IC and ID) for a test run with initial translator 
and gimbal position offsets. In the case of the translator, 
the position offsets are analogous to and indistinguish- 
able from center-of-mass misalignments in the x- and 
(/-direction. Initial values were y T — 0.0213 m (0.07 ft), 
x T = —0.0244 m (—0.08 ft), and a = 0.095 rad. 

The transient responses of the x- and (/-axes, as con- 
trolled by the (/-translator and x-translator, respectively, 
were nicely damped and therefore quite satisfactory. The 
“ghtch” in y T ( t ) at about 310 s is as yet unexplained, 
although there is some suspicion that it was the result of 
a power switching transient to the ion engine bank operat- 
ing in the vacuum chamber. (Such power transients had 
been observed once or twice earlier to couple enough 
energy to some TVCS cables to momentarily send a short 
burst of pulses to one or another TVC actuator.) Note also 
the roughly 10 times greater y-axis response magnitudes 
than in the x-axis, a consequence, as mentioned earlier, of 
dropping the (/-axis sensor gain by a factor of 10 to com- 


pensate for about 10 times less inertia. VCO and compen- 
sation loop parameters were the same for both translator 
loops. 


The gimbaled z-axis response obviously was not well 
damped and therefore somewhat less than satisfactory in 
performance. A look at the linearized z-axis loop root loci 
in Fig. 16 shows the actual operating point for the test run 
results of Figs. 14 and 15. With either one or both outer 
gimbaled engines operating, the damping factor on the 
main locus is a rather poor 0.24. The root loci shown cor- 
respond to the gimbal loop control parameters used in the 
test simulation, which are: 


t x = 500 s 
t 2 = 440 s 
Tg = 0.1 s 
Kx = l.OV/step/s 


K 2 = 3.65 
K F = 0.50 
K vc0 — 50 step/s/V 
K s = 295 V/rad 


(The translator loop root loci are as shown in Fig. 4.) It is 
obvious that the value of K 2 in particular was poorly 
chosen for the gimbal loop, and K? could also be modi- 
fied somewhat to good effect. In fact, by using K t = 7.0, 
K F = 0.35, and all other parameters unchanged, a much 
improved operating point in the S-plane is produced, 
as shown in Fig. 17. [For the gimbal loop, F = 0.057 
N-m/rad (0.0424 ft Ib/rad) and K» = 0.97 (10 *) rad/ 
step.] The transient response for such a modification will 
be shown in the Section IH digital computer simulations. 
However, no adjustments of that kind were made in the 
TVCS gimbal loop electronic hardware, and the remain- 
ing test run results reflect the poorly damped gimbal 
responses. 


Figure 18 is another recording of the TVCS response to 
initial translator and gimbal displacements, showing the 
pitch phase-plane transient. Of particular interest in this 
x-y plot is the very evident effect of analog computer inte- 
grator drift during the test. When the phase-plane tra- 
jectory (which starts at 9 X = 0, 8 fx = 0) first crosses the 
9 X = 0 line at about t = 470 s, the requirement that the 
slope of the curve be infinite is very nearly met, although 
there is a slight positive bias (drift) in the rate. However, 
at the next apparent crossing ( t ss 1500 s), the proper rate/ 
position relationships are obviously not present, an indica- 
tion of significant levels of drift in both the rate and posi- 
tion integrators. 

One further clarification of the results shown in Figs. 14, 
15, and 18 is required. Vehicle bus rotations about the x- 
and z-axes are labeled 0 /x and 8 /x , respectively (see Fig. 13) 
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and do include the effect of solar array vibrations, although 
no effects are visible. However, because of the manner in 
which the solar array transfer functions happened to be 
programmed (see Fig. 13), i.e., inserted between position 
variables 8 X and 8 fx (or 8 Z and $ /2 ), the true bus angular 
rates 0 /x and 6 fz were not available. As a result, 8 X and 9 Z 
were recorded instead, actually representing the bus’ rota- 
tional rates minus the perturbing vibrations of the array 
appendages. These perturbations in rate, though, were 
quite small and will be shown in the digital computer 
simulations of Section IH. The variables 8 y and 0 V are 
strictly rigid-body responses, since appendage modes of 
vibration were not included in this axis. 

At the conclusion of the TVCS transient test response to 
initial engine displacements (t = 2500 s), an engine-out 
condition was simulated by instantaneously setting f 2 = 0. 
The TVCS response to that action was recorded, and the 
actual strip charts are shown in Figs. 19 and 20. 

Engine gimbal angle and roll angular position remained 
reasonably quiescent during the engine-out transient, as 
they should, in spite of the 50% reduction in effective roll 
control torque. Translator displacements [each is required 
to move 0.114 m (0.375 ft) to put the resultant thrust vec- 
tor through the vehicle c.m.] were well damped, as pre- 
dicted by the linear analysis (Fig. 3). And response time 
of the translator loops was predictably slowed as a result 
of the smaller control forces and the accompanying lower 
loop gains. Total translator excursions during the transient 
were less than 0.183 m (0.6 ft) out of a total available 
travel of 0.67 m (2.2 ft). The peak pitch angle excursion of 
the vehicle during the transient was 2.3 deg. 

G. Test Problems 

Because of the likelihood of continued TVCS hardware 
tests occurring in the immediate future, it may prove valu- 
able to recount some of the operational difficulties expe- 
rienced thus far in the hope that at lea£t some of these 
may be overcome. 

Insofar as the TVCS portions of the SEPST III tests 
are concerned, the simulations were strictly a one-man, 
part-time (approximately half-time) operation spread over 
roughly a 4-month period. At the same time, a great deal 
of effort in terms of man-hours was understandably de- 
voted to ion engine vacuum chamber tests with relatively 
frequent test shutdowns and system modifications. As a 
consequence, TVCS closed-loop -simulations were, at 
times, precluded because engine mechanisms could not 
be moved while the modifications were being effected. 


While these types of conflicts are inevitable, they do sig- 
nificantly reduce the amount of useful test time and should 
be planned for. 

Without a doubt the most frustrating part of the TVCS 
tests was simply a consequence of the system’s own char- 
acteristics. This was the extremely long period of time 
(usually at least 40 min) necessary to generate just a single 
transient response from the system. Of course, hand in 
hand with the slow, real-time TVCS response came the 
problem of analog integrator drift, which, after 10-15 
min of computation time, began to destroy any accu- 
rate position-rate relationship in the vehicle dynamics 
model. The drift degradation of the simulations made 
them useless for predicting steady-state angular rates or 
positions of the vehicle. It also goes without saying that 
the bandwidth of this real-time model made any attempts 
to use it repetitively as an analytical tool impossible. 

The time factor difficulties mentioned above were mag- 
nified by a variety of other problems, some of which were 
probably to be expected. One of these was related pri- 
marily to the age of the EAI 231-R analog computer. Now 
12 years old, the computer exhibits a number of ills, not 
the least of which are intermittent patch-board connec- 
tions. Oxidized, dirty, and sometimes bent contacts and 
patch-cord pins are the culprits, and it requires some 
patient searching and experimentation to find the trouble 
and achieve a reliable patch-board setup. 

A faulty electrical connection at the signal patch-board 
was also the source of lost test time. Several test runs 
were interrupted by either intermittent or total dropout 
of the x-translator output signal to the analog computer. 
Although at first the location of the break was unknown, 
it was ultimately traced to a bad solder connection behind 
the board. Because of its location, however, the connec- 
tion was not repairable during the active SEPST III 
engine test period, and it was necessary to live with the 
problem for TVCS testing. 

Since the simulation control and data recording were 
done at a site remote from the location of the TVCS hard- 
ware, particularly the electronics panel, it was extremely 
inconvenient and time-consuming to travel back and forth 
to accomplish such simple tasks as turning the control 
electronics on and off, resetting or disabling the feedback 
compensation circuits, or initializing the positions^of trans- 
lator or gimbal mechanisms. It would be quite valuable in 
future testing to have some capability for remote and auto- 
matic control of the system’s configuration and initial state. 
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In general, efficient future testing, at least from the 
TVCS standpoint, would appear to require better coordi- 
nation of diverse and often conflicting test objectives and 
a greatly enhanced capability for automatically controlling 
the test equipment configuration and operation. 

H. Digital Computer Simulations of the Solar 
Electric TVCS 

As a check on the SEPST III TVCS test results, particu- 
larly those shown in Figs. 14, 15, 18, 19, and 20, the system 
equations were programmed for digital computer solution 
using Continuous System Simulation Language (CSSL III) 
on the Univac 1108. The vehicle dynamic characteristics 
were programmed exactly as they were for the analog 
computer and as they have been described in Sections IA- 
IE. Control parameter values were, for the most part, 
taken at their nominal values (e.g., actuator gear ratios, 
sensor gains) except for a few compensation loop param- 
eter changes of the order of 5% that were made to better 
reproduce actual test responses. 

The digital simulation included only one significant 
nonlinearity (except for limits on translator/gimbal 
excursions-limits that were never reached) not shown in 
the typical control loop block diagram of Fig. 1. This 
involved a somewhat more accurate simulation of the 
actual circuit implementation of the time constant t, and 
associated gain. 

The circuit of interest is shown in Fig. 12d, where t x is 
mechanized using the 100-/xf capacitor at the input of 
operational amplifier A. At very low pulse rates (low duty 
cycle), the impedance of the circuit from the capacitor 
back through the transistor switch is primarily that of a 
tumed-off switch, or, in other words, infinite. Therefore, 
the effective time constant is that of the 5-MO input resis- 
tor and 100-/J capacitor, or 500 s (a number which is used 
for r, in all the root locus calculations shown in this re- 
port). At the other extreme, at the highest pulse rates of 
the VCO (100% duty cycle), the impedance is that of a 
turned-on switch in series with 5 Mq, so that = 250 s. 
(Of course, placing the capacitor in the feedback of ampli- 
fier A would have prevented these complications, but the 
problem was not recognized until tests were completed.) 
Thus, strictly speaking, r t is a function of VCO pulse rate 
(although n = 500 s is a good approximation) and so is 
the dc gain through amplifier A. Figure 21 shows the 
actual functional dependence in the desired transfer func- 
tion equivalent assuming 25-V pulses, 0.20 s wide, at the 
transistor switch output. A straight-line approximation to 
the function shown in Fig. 21, very accurate over the pulse 


frequency range 0-10 pulses/s, was included in the digi- 
tal simulation. 

Figures 22-32 show the results of a digital simulation 
of the TVCS test run for the initial engine offsets described 
in Section IF and in Figs. 14, 15, and 18. Translator and 
gimbal control loop parameters are those given previously 
and are consistent with the root loci in Figs. 3 and 16. 
A program listing for the digital simulation is given in 
Fig. 33. All transient peaks and their times of occurrence 
agree quite well with the TVCS test recordings. However, 
one noticeable difference occurs in the rate of angular 
position (0/1,01/,0/z) error reduction over the 40+ min 
simulation period. The analog computer solutions (Figs. 14 
and 15) show a much slower reduction, which, as indicated 
before, is very probably due in large part to the cumulative 
effect of drift in both the rate and position integrators. 
There may also be some nonlinear effects occurring in the 
near-steady-state, limit-cycle-type operation of the control 
loops, particularly in the VCO (such as bias voltages at the 
input) feedback compensation loop, and actuators which 
have not been modeled in enough detail. Remember that 
the digital program used does not actually deal with dis- 
crete pulses but idealizes all system variables as smooth, 
continuous functions. While the former level of detail 
could be added, the digital simulation would run at least 
10 times more slowly (and 10 times more costly) as a result 
of its efforts to correctly integrate the effects of thousands 
of narrow VCO pulses and stepper motor steps. 

Notice also the plot of 6 tx (t) in Fig. 23 as compared to 
6 X in Fig. 22. Except for some very small solar array vibra- 
tion during the first 4 min, the responses are identical. The 
plot of 6f Z (Fig. 29) shows no noticeable vibration at all. 

In Section IF, it was mentioned that a substantial im- 
provement in gimbal loop response could be made by 
using the values fC 2 = 7.0 and K F = 0.35. Figures 34-36, 
results of a digital simulation, demonstrate the improve- 
ment over the z-axis responses of Figs. 15 and 29-31. 
Translator loop transient responses were unchanged, since 
the control axes are effectively independent. Note, how- 
ever, that the response curve for 6/~ now shows some per- 
ceptible solar array vibration. 

A digital computer simulation was also performed for 
the “engine-out” TVCS test described in Section IF (see 
Figs. 19 and 20). Plots of the computer results are given 
in Figs. 37-45. This simulation is actually a continuation 
of the TVCS response to initial translator and gimbal dis- 
placements with f 2 set equal to zero at t = 2450 s, just as 
was done in the actual SEPST III test run. Again, a com- 
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parison of the digital computer solutions to the actual test 
run recordings of Figs. 19 and 20 will show very good 
agreement throughout the transient period. 

I. A Digital Controller and Vehicle Dynamics Simulator 
for SEPST Testing 

Throughout the SEPST III tests, a PDP-11 digital com- 
puter was used as a test controller (see Refs. 9 and 10), 
performing such functions as turning engine power- 
conditioning equipment on and off, monitoring engine 
thrust performance, engine failure detection, adjusting 
thrust levels in relation to pre-stored “flight plans,” and 
simulating various ground commands. Virtually no control 
was exercised, however, over the functioning of the TVCS 
electronics, actuators, or analog computer simulation of 
vehicle dynamics. But the PDP-11 does appear to have the 
computational capability not only to interact with the 
TVCS electronics (i.e., for gain switching or monitoring 
purposes) but actually to replace the electronics and per- 
form the thrust vector control function itself. 

Space has already been provided in the PDP-ll’s 8K 
memory for a prototype subprogram which has been writ- 
ten to perform the per-axis control function shown in 
Fig. 46. Just like the present TVCS electronics, the com- 
puter will receive a simulated celestial sensor error voltage 
and will derive a sequence of actuator drive pulses as a 
function of both the sensor error and a feedback com- 
pensation signal. 

The computer will be capable of updating the pulse 
output to the actuator drivers every 20 ms. A very simple 
integration algorithm may be used (such as Euler’s Method 
shown in Fig. 46), since the iteration rate will be extremely 
fast compared to system time constants. 

A basic algorithm for duplicating the TVCS electronics 
analog functions on a digital computer is shown in Fig. 47. 
Initial attempts at performing the control function with 
the PDP-11 will be limited to those functions shown in 
Fig. 46, first for a single axis and then, hopefully, for all 
three axes. However, if these tests prove successful and 
computer storage space permits, the PDP-ll’s task will be 
enlarged to include the simulation of vehicle rigid-body 
dynamics as well (two additional integrations per con- 
trol axis), including the detailed engine thrust-torque 
equations. 

If the computer is able to handle these functions, a num- 
ber of current testing problems would be effectively 


solved. Certainly the problem of analog integrator drift 
during the long TVCS test runs would be alleviated. 
Absolute control could easily be maintained over almost 
all initial conditions of the system variables, and even the 
initial translator and gimbal positions could be set by com- 
puter command rather simply. Complete flexibility with 
respect to system parameters would be provided; changes 
would be made by a simple typewriter entry. Even the 
control law(s) would be changed relatively quickly with 
only a software modification. Further, the concept of an 
on-board programmable processor for use as a spacecraft 
attitude controller has already received much support, 
particularly for extended-life missions (e.g., to the outer 
planets). One design, the Hybrid Programmable Attitude 
Control Electronics (HYPACE), is already in the develop- 
mental testing stage (Ref. 9) and could be used, among 
other things, to try out promising solar electric thrust 
vector control law software with real actuators. 


II. A Solar Electric TVCS Using Rate Estimation 

A. Linearized System Model 

The utility of a state estimation circuit, derived by Kopf 
(Ref. 12) specifically for the task of estimating the angular 
rotation rate of a space vehicle, was investigated as it 
might apply to the same basic solar electric thrust vector 
control problem discussed in Section I. The estimator, 
receiving celestial sensor measurements of the space 
vehicle’s angular position, both filters the sensor signal 
and derives an estimate of vehicle angular rate. 

A transfer function block diagram of the proposed 
single-axis translator (or gimbal) control loop, with esti- 
mator, is shown in Fig. 48. Critical parameters associated 
with the estimator are the filter “time constant” B and the 
rate-to-position gain K R . To the extent that the estimator 
represents a model of the vehicle rigid-body dynamics with 
corrective feedback through the gains 1/B and 1/2B 2 , one 
has the option of improving the estimation accuracy by 
introducing a measure of the known applied control 
torques on the vehicle T c , divided by the proper inertia 
constant. However, in the analyses and simulations that 
follow, this option is foregone as an unnecessary compli- 
cation. If one sums the rate and position estimates into a 
single output, the estimator’s transfer function has the 
following form: 


V„ (S) S 2 (2 B 2 + BK r ) + S(K r + 3 B) + 1 
V s (S) (BS + 1) (2B 2 S 2 + 2BS + 1) 
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Factoring the numerator polynomial, the roots are 

K R +2B 1 

ai ~ 2B 2 + BK r ’ 2 2B + K R 

Assuming that K s > > B, then 



and 

Vo(S)_ (K r S + 1) 

V s (S) ~2 B 2 S 2 + 2BS + 1 

If the actual sensor position signal, rather than the posi- 
tion estimate, is summed with the rate estimate, the trans- 
fer function becomes 

Vo (S) 2 B 2 S 2 + S (K R + 2 B) + 1 

Vs (S) ~ 2B 2 S 2 + 2BS + 1 


If K b > > B, the numerator roots are given approxi- 
mately by 



and 

/2B 2 \ 

y.m < K » s + 1 >Ur s + 1 ) 

V g (S) = 2B 2 S 2 + 2BS + 1 

In this application, the filtered position signal V s will be 
used to aid in minimizing any potential noise problems at 
the VCO input. The linearized single-axis open-loop trans- 
fer function for the controller with rate estimator (Fig. 48) 
becomes 


T c (S) _ FKuKg K v co (r,S + 1) (K*S + 1) (S 2 + 2|<oS + « 2 ) 

T (S) IS 3 (t»S + 1) ( ri S + KpKvco + I) (2 B 2 S 2 + 2BS + 1) (RS 2 + 2^S + o, 2 ) 


The closed-loop root loci for the single-axis controller 
shown in Fig. 48 were calculated for a number of param- 
eter variations, and a sample of the results is given in 
Fig. 49. The parameters used to calculate the loci in Fig. 49 
were as follows: 

B = 1 
Kr = 200 s 
Kvco = 50pulses/s/V 
Kr = 0.20 
n = 200 s 
Ks = 295 V/rad 

K m = 0.64 X 10' 4 m/step (0.00021 ft/step) 

F = 0.089-0.267 N (0.02-0.06 lb) 

I = 20,337 kg-m 2 (15,000 slug-ft 2 ) 

The celestial sensor time constant was ignored, as were 
the poles and zeros contributed by the fundamental solar 


array vibration mode. Good damping is obtained over the 
entire range of engine thrusts, which was not the case 
when the controller described in Section I was used. 
Transient response of the loop, using the translator param- 
eters listed above, to a 0.061-m (0.2-ft) vehicle c.m. offset 
is shown in Figs. 56-52. The speed of response with the 
estimator is at least twice that shown in the responses, say 
in Figs. 23-25, which are for the present controller with 
less than half the disturbance. Also, the peak angular 
position transient with the estimator is 1.35 mrad for a 
0.061-m (0.2-ft) c.m. offset. With the present controller and 
a 0.0213-m (0.07-ft) offset, 0 error peaks at over 7 mrad. 


B. Three-Axis System Simulations 

In order to get a more direct comparison between the 
rate estimating controller and the present scheme, the 
digital simulation of the complete three-axis system with 
solar array vibration, used in Section IH, was employed. 
VCO gain in all three axes was identical at K vc o = 
50 pulses /s/V, as was the estimator parameter, B — 1. 
However, the control parameters K fi , r 1; and K F were 
tailored to provide good (although not really optimal) 
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responses in each control axis. The values used are given 
in Table 3. 


Table 3. Rate estimator control parameters for 
three-axis SEMMS vehicle simulation 


Parameter 

x-axis 

y-axis 

z-axis 

Kr, s 

200 

100 

250 

T lt S 

200 

100 

250 

Kf, V/pulse/s 

0.20 

0.30 

0.175 

Ks, V/rad 

250 

250 

295 


Figures 53-61 show system responses to the identical 
initial conditions present in the SEPST III test and cor- 
responding digital simulation, i.e., initial translator and 
gimbal position offsets. Comparisons of x- and z-axis re- 
sponses in particular show the simulated rate estimator 
loops to be well damped and substantially faster than the 
present control hardware. In fact, the widened bandwidth 
in the estimator loop is especially evident in the x-axis 
response, where solar array perturbations seem to be in 
somewhat more evidence than before. The gimbal loop 
(z-axis) response has been drastically improved over that 
given in Section I as a result of the estimator and proper 
control parameter choices. Notice also that the low-inertia, 
yaw axis control loop has had its control parameters 
( K r , ti, Kf) adjusted to provide good response without 
lowering the celestial sensor gain K s , as was done in 
Section I. 

To complete the comparisons, the rate-estimating TVCS 
was also simulated for the same engine-out condition 
described in Section I. The resulting transient responses 
are given in Figs. 62-70. In pitch and yaw (x and y), the 


translator-controlled axes, recovery from the single-engine 
shut-down (f 2 = 0) is fast and nicely damped. The gimbal 
axis transient is not nearly as well damped but extremely 
small in magnitude. (The rate-estimating TVCS simulation 
program is listed in Fig. 71.) 

In summary, the rate-estimating solar electric TVCS 
certainly promises a vastly improved control capability 
over the present approach. A control bandwidth improve- 
ment by a factor of 2-3 does not seem unreasonable, and 
yet this still does not threaten to result in any significant 
solar array interaction problems, much less any serious 
destabilizing effects. 

C. Future Work 

The immediate future will see efforts to bench-test a 
single-axis version of the rate-estimator control loop using 
real actuators and simulated rigid-vehicle attitude dynam- 
ics. It will be necessary to demonstrate that rate-to-position 
gains on the order of 200-250 s are practical to achieve 
with conventional circuitry so that these gains, in combina- 
tion with an estimator time constant of about 1.0 s, do not 
present noise problems or excessive actuator limit-cycle 
stepping rates. It does appear, however, that the value of 
B could be increased substantially if necessary to obtain 
greater filtering action without significantly altering the 
system transient response. 

In addition, an attempt will be made to program an 
estimator-TVCS loop on the PDP-11 computer to see if an 
analog circuit implementation can be dispensed with 
entirely, along the same lines and for the same reasons 
discussed in Section I. 
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Fig. 3. Closed-loop root loci for the single-axis TVCS 
(rigid vehicle) 
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Fig. 4. SEMMS vehicle three-engine configuration 
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Fig. 5. SEPST III TVCS test configuration: (a) thruster array, (b) gimbal actuator, (c) translator actuator, 
(d) closed-loop control panel, (e) thrust vector control electronics breadboard 
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Fig. 12a. TVCS feedback compensation and VCO output circuitry 
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Fig. 12b. Translator control electronics 
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Fig. 12c. Gimbal control electronics 



Fig. 12d. Modified TVCS feedback compensation and 
VCO output circuitry 
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Fig. 13. SEPST III TVCS test analog computer program 
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Fig. 16. Gimbaled axis closed-loop root loci for the TVCS test 


Fig. 17. Gimbaled axis closed-loop root loci using modified 
control parameters 
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Fig. 18. SEPST III spacecraft pitch axis response to 
0.02-m (0.07-ft) CG offset 
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Fig. 19. TVCS test responses to engine-out (f 2 = 0) transient condition (e x , e/i, yr, e y ) 
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Fig. 20. TVCS test responses to engine-out (f 2 = 0) transient condition (e v , xr, 6/z, a) 
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Fig. 25. yr vs. time, computer-simulated test response to 
initial translator and gimbal displacements 
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Fig. 26. by vs. time, computer-simulated test response to 
initial translator and gimbal displacements 
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Fig. 27. e y vs. time, computer-simulated test response to 
initial translator and gimbal displacements 



Fig. 28. xt vs. time, computer-simulated test response to 
initial translator and gimbal displacements 
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Fig. 31. a vs. time, computer-simulated test response to 
initial translator and gimbal displacements 


Fig. 29. hu vs. time, computer-simulated test response to 
initial translator and gimbal displacements 
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Fig. 30. 6f Z vs. time, computer-simulated test response to 
initial translator and gimbal displacements 
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Fig. 32. e z vs. 8 ix, computer-simulated test response to 
initial translator and gimbal displacements 
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Fig. 53. $tx vs. time, simulated TVCS (with estimator) response 
to initial engine array offsets 
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Fig. 54. Six vs. time, simulated TVCS (with estimator) response 
to initial engine array offsets 
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Fig. 55. yr vs. time, simulated TVCS (with estimator) response 
to initial engine array offsets 
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Fig. 56. 6 y vs. time, simulated TVCS (with estimator) response 
to initial engine array offsets 
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Fig. 57. 9 y vs. time, simulated TVCS (with estimator) response 
to initial engine array offsets 
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Fig. 58. xt vs. time, simulated TVCS (with estimator) response 
to initial engine array offsets 
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Fig. 59. o f z vs. time, simulated TVCS (with estimator) response Fig. 60. $iz vs. time, simulated TVCS (with estimator) response 
to initial engine array offsets to initial engine array offsets 
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Fig. 61. a vs. time, simulated TVCS (with estimator) response 
to initial engine array offsets 
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Fig. 62. 6j Z vs. time, simulated TVCS (with estimator) response Fig. 63. etz vs. time, simulated TVCS (with estimator) response 
to engine-out transient to engine-out transient 
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Fig. 64. yr vs. time, simulated TVCS (with estimator) response 
to engine-out transient 
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Fig. 65. i v vs. time, simulated TVCS (with estimator) response Fi 8- 66 vs - time - Emulated TVCS (with estimator) response 
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Fig. 67. xt vs. time, simulated TVCS (with estimator) response 
to engine-out transient 
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Fig. 68. S/z vs. time, simulated TVCS (with estimator) response 
to engine-out transient 
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Fig. 69. 6/t vs. time, simulated TVCS (with estimator) response 
to engine-out transient 
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Fig. 70. a vs. time, simulated TVCS (with estimator) response 

_ to engine-out transient - 
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program solar electric tvc (with rate estimator) 

TX= CF1+F2+F3 ) * < YT-YCG ) + ( F3-F2 ) * < A+ALP*D*DEL ) 
TY=(Fl+F2+F3)*(xCG-XT)+(F3-F2>* ( A“ALP*D*DEL ) 
TZ=(F3-F2 >*<XT-XCG-YT+YCG)*ALP*DeL-(F 2+F3)*2.*A*ALP*DEL 
THXDD=< Ill*TX+I12*TY+I13*T2)*1000. 

THYDD=< Il2*TX+l22*TY+I23*TZ)*1000. 

THZDD=Ul3*TX+I23*TY + I33*TZ)*1000. 

THXO= INTEG ( THXDO t THXDO ) 

THY0=INT£G(THY0D*THY00) 

THZD= INTE6 ( THZDD » THZDO ) 

thy=integ<thyd»o.) 

XXDD= ( THXD-XX* (S1GX**2) -2 . *ZETX*SlGX*XXD ) /RX 
XZ0D=(THZD-XZ*(SIGZ**2)-2.*ZETZ*SIGZ*XZD)/RZ 
XXD=1NTEG<XXDD»0.) 

XZD=INTEG(XZDOf U. ) 

XX=lNTEG(XXD»0. ) 

XZ=lNTEG(XZO»0.) 

ThFXD=XX0D+2.*ZETX*SIGX*XXD+(SIGX»*2) *xx 
ThFZ0=XZD0+2.*ZETZ*SI6Z*XZD+(SIGZ**2) *XZ 
THFX=INTEG(THFXD*0. ) 

THFZ=INTEG(THFZD»0. ) 

THFXOE=InTEG( <THFX-THFXE)/(2.*<B**2> > *0. ) 

THFXE=INTEG( (THFXOE-(THFXE-THFX) /B) (0. ) 

THFZDE=InTEG< (THFZ-THFZE)/(2.*(8**2) > ► 0 . ) 

THFZE=INTEG( (THFZDE-<THFZE-THFZ)/B) *0. ) 

THYDE=INTEG( (THY-THYE)/(2.*(B**2> > *0. ) 

THYE=INTEG( (THYLE-(THYE-THY)/B) »o. ) 

VXE=KSX* ( THFXE+KRX+THFXOE ) 

VYE=KSY*(THYE+KRY*THYDE) 

VZE=KSZ*(THFZE+KRZ*THFZDE) 

VC0XT = B0UND (-50. » 50 • *KVCOXT* (VYE-VFXT ) ) 

VCOYT=BOUND(-50. *50. * KVCOYT* ( VXE-VFYT > ) 

VCOG=BOUND ( -50 . , 50 . * K VCOG* ( VZE-VFG ) ) 

VFXT=INTEG( (KFXT*VCOXT-VFXT)/T1XT*0. ) 

VFYT=INT£G( (KFYT*VCOYT-VFYT)/T1YT>0. ) 

VFG=INTEG( (KFG*\iCOG-VFG)/T1G*0. ) 

XT=LlMlNT(KMXT*VCOXTfXTOf-l.l*l.l) 

YT=LIMlNT(-KMYT*VCOYT*YTO*-l.l*l.l) 

ALP=LIMINT(KMA*VC0G* ALPOf-,2* .2) 

CINTERVAL CI=5.0 
STPCLK 200. 

TERMT (T.GE.FlNTM) 

CONSTANT B=l. *KFXT=.3U*KFYT=.2U*KFG=.175 
CONSTANT T1 XT=i00. »T1YT=200. »T1G=250. 

CONSTANT KRX=200. »KRY=100. *KRZ=25Q. 
constant fi=.020»f3=.020*f2=.020 
CONSTANT xcg=o. »ycg=o. 
constant a=.75»d=4. *bel=.707i 

CONSTANT IU=.709089E-4* I12=-.95843E-7r U3=-.292855E-7 
CONSTANT I22=.658164E-3»I23=-.202686E-5» 1 33= . 785683E-4 
CONSTANT THXOO=0. «THYD0=0. t THZDO=0. 

CONSTANT KSX=.250rKSY=.2500»KSZ=.295 
CONSTANT KVCOXT=50, ?KVCOYT=5.0.. *KVC0G=5fl. 

CONSTANT KMXT=2 • 06E-4 » KMYT=2 . 08E-4 » KMA= . 97E-4 
CONSTANT XTO=-. 0800 *YTO=. 070 »ALPO=. 095 
CONSTANT RX=.0983»RZ=. 03413 
CONSTANT SIGX=.3027*SiGZ=.2487 
CONSTANT ZETX=.005*ZETZ=.005 
CONSTANT FINTM=2450. 

OUTPUT 20 . XXO » THY , XZO » THXD . THYo r THZD * XT » YT , ALP > vCOxT * VCoYT > ... 
ThFXD»THFZD*VCOG» THFX.THFZ 

PREPAR XT* YT *ALP»THFXO#THYD*THFZ»ThFX»THY f THFZD'THXO 

END 


Fig. 71. Simulation program for rate-estimating TVCS 
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